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Abstract 

Kaon condensation may play an important role in the structure of hadronic matter at densities 
greater than that of nuclear matter, as exist in the interior of neutron stars. We present the 
results of the first lattice QCD investigation of kaon condensation obtained by studying systems 
containing up to twelve negatively charged kaons. Surprisingly, the properties of the condensate 
that we calculate are remarkably well reproduced by leading order chiral perturbation theory. In 
our analysis, we also determine the three-kaon interaction from the multi-kaon systems and update 
our results for pion condensates. 



I. INTRODUCTION 

A key ingredient in determining whether supernovae evolve into black-holes or neutron stars 
is the nuclear equation of state (NEOS) [1J. The NEOS depends upon the degrees of freedom 
at a given density, which in turn are determined by the masses of the various hadrons in the 
medium and their interactions. One possibility is that it is more energetically favorable for 
a Bose-Einstein condensate of negatively charged mesons to form in the hadronic matter. 
Naively, one expects pion condensation to occur before kaon condensation 1 , however, the 
repulsive pion-nucleon interactions and the attractive kaon-nucleon interaction lead to kaon 
condensation being the prime candidate for softening the NEOS at moderate nuclear den- 
sities. In theoretical explorations of kaon condensation, chiral perturbation theory (%PT) 
is used to describe the kaon-nucleon interactions and the kaon self-interactions at leading 
order (LO) or next-to-leading order (NLO) in the chiral expansion [3]. The energy of the 
hadronic plus electronic system is minimized with respect to the vacuum expectation values 
(the condensates) of the meson fields subject to the constraints imposed by charge neutral- 
ity and the conservation of baryon number. In the case of non-interacting mesons, once 
the chemical potential exceeds the rest mass of a given species of meson a condensate will 
form with an infinite density of mesons. However, in the case of a meson with repulsive 
self-interactions, the number density of mesons will remain finite and a smooth function of 
the chemical potential. Depending upon the many-body technique employed to describe the 
system, a condensate of K~ ] s may become energetically favorable at baryon densities, n b , 
as low as three times nuclear matter density (for a nice discussion, see ref. [2]). 

In this work we perform the first Lattice QCD studies of kaon condensation. In particular, 
we calculate the correlation functions of systems containing up to twelve (n = 12) K~'s on 
the coarse MILC lattices with spatial extent L ~ 2.5 fm and L ~ 3.5 fm, both with lattice 
spacing b ~ 0.125 fm, and on the fine MILC lattices with spatial extent L ~ 2.5 fm and 
lattice spacing b ~ 0.090 fm. LO %PT is found to describe the chemical potential, and its 
density dependence, remarkably well over the densities explored. This work is an extension 
of our previous work on charged-pion condensates [5J E] and, for completeness, we update 
that analysis with the new ensembles studied herein. 



II. CONDENSATION OF PSEUDO-GOLDSTONE BOSONS IN X PT 

At LO in xPT> the Lagrange density describing the low-energy dynamics of the pseudo- 
Goldstone bosons associated with the spontaneous breaking of the approximate SU(3)l <2> 
SU(3)r chiral symmetry of QCD has the form 

C = 4 Tr [ <9 M E d^ ] + A Tr [ m q (E + E 1 " ) ] , (1) 



i 



A second possibility is that it is energetically favorable to have S~'s present in the nuclear material in 
addition to neutrons and protons. The in-medium mass of the £~ is expected to be significantly less 
than in free-space due to its interactions with neutrons, but the precise mass-shift is uncertain due to the 
model dependence of existing theoretical calculations |2J. Experimentally, little is known about the S _ - 
neutron interaction and a precise Lattice QCD calculation of these interactions [5] would greatly reduce 
the uncertainty in this particular contribution to the NEOS. 



where £ is the exponential of the meson fields, 
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and where / is the chiral limit of the pion decay constant (in this normalization convention 
the physical pion decay constant is f n ~ 132 MeV). At LO in the chiral expansion, the 
quark mass matrix, m q , can be written in terms of the meson masses, 
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The Lagrange density in eq. (fij), in addition to describing the free-field dynamics of the 
mesons, also contains interactions containing arbitrary numbers of mesons. At zero momen- 
tum, the LO contributions to the interactions between three charged pions or between three 
charged kaons are given by 
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This contribution provides the naive dimensional analysis (NDA) estimate of the size of the 
three-meson interaction [3, E] of l/(m 7r /' 4 ) for the tc~7t~7t~ interaction and l/(rriKf 4 ) for 
the K~K~K~ interaction in the non-relativistic theory describing the very- low momentum 
interactions between the mesons 2 . 

In the presence of a chemical potential for the third component of isospin, ///, and 
strangeness, fj,$, the partial derivatives in eq. ([I]) become covariant derivatives [U QUI [TTJ 
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In the situation where //„.- = — /xj > m n and jj,k- = /As — \ni < rn K , it is energetically 
favorable to form a n~ condensate but not a K~ condensate, and at LO in xPT the number 
density of vr-'s is HU] 
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The analysis leading to this result easily extended to kaon-condensates [TTJ. For //„.- < m n 
and fix- > m K there will be a condensate of K~'s but not of 7r~'s 3 . For chemical potentials 



2 The relativistic and non-relativistic fields differ in their normalization by a factor of y/rn. 

3 Mixed states with different numbers of pions and kaons are the subject of ongoing work. 



in this range, the number density of K 's as a function of the chemical potential, fix- 
at LO in ;\PT, 
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In eqs. ([7]) and (J8J), we have used LO x?T relations to express the result in terms of the 
physical masses and decay constants. It is this relation that we will compare with Lattice 
QCD calculations. The reason for "testing" the relation in eq. (J8J) is that, given the relatively 
large expansion parameters of SU(3)l <8> SU(3)r x?T, it is natural to expect that there 
will be corrections to Pk-, and hence to the dependence of jjlk- on px-, at the level of 
m 2 K / 'A 2 ~ 25%, where A x ~ 1 GeV is the scale of chiral symmetry breaking. 

In order to assess the relevance of the present work, it is instructive to compare the 
value of the ^"-condensates we calculate to those found to be relevant to the dense nuclear 
matter in the interior of neutron stars. Decomposing the charged kaon fields into their 
real and imaginary parts, K ± = (K\ =F 1K2) /v2, and denoting the condensates (Ki) = V\ 
and (K2) = V2 and further defining (K) = \Jv\ + 7jf /\/2, the condensate is related to the 
chemical potential via 
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Condensates of magnitude (K)/ f K < 0.3 are found to be possible at the densities in the 
interior of neutron stars [I]. This corresponds to //#- /thk — 1 ^ 0.1, which is contained in 
the region explored in the Lattice QCD calculations described in this work. 

III. MULTI-MESON ENERGIES IN A FINITE VOLUME 

In recent works pj |H1 H2], the analytic volume dependence of the ground-state energy of n 
identical bosons in a periodic volume of side L has been computed to 0(L~ 7 ), extending 
the classic results of Bogoliubov [13] and Lee, Huang and Yang [14J . The resulting shift in 
ground-state energy of n particles of mass M due to their interactions is [8] 
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where the parameter a is related to the scattering length , a, and the effective range, r, by 

2 > (i-(±)i) ■ (id 
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The geometric constants that enter into eq. (10) are 
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(12) 
and n C m = n\/m\/(n — m)\. The three-body contribution to the energy-shift given in 



eq. (10) is represented by the parameter rj 3 , which is a combination of the volume-dependent, 
renormalization group invariant quantity, r\\ , and contributions from the two-body scattering 
length and effective range, 

-l _ L ( , „ ( a \ ^\ 727ra 4 r 



S = 'H 1 - 6 UJ x J + ^r I • < 13 > 

where 

r , . 647ra 4 ( r- \ , . ,, 96a 4 _ . , 

t = V3(ri + - jr (sV3-An) log(/iL) - ^S MS ■ (14) 

The quantity 773 (/V) is the renormalization scale dependent coefficient of the three- ir + in- 
teraction that appears in the effective Hamiltonian density describing the system [Sj. The 
quantity S is renormalization scheme dependent and its value in the minimal subtraction 
(MS) scheme is S MS = -185.12506. 



For n = 2, the last two terms in eq. (10) vanish and the remaining terms constitute the 



small a/L expansion of the exact eigenvalue equation derived by Luscher [151 HH] , 

pcot6(p) = isf f^V ] , (15) 



which is valid below the inelastic threshold, where p cot 5(p) is the real part of the inverse 
scattering amplitude. The regulated three-dimensional sum is [151 EH E] 

|J|<A 

S(x) = J2n~2 47rA > ( 16 ) 

. |J| X 

where the summation is over all triplets of integers j such that | j | < A and the limit A — > 00 
is implicit. 



4 In this work we use the Nuclear Physics sign convention for the scattering length, opposite to that used in 
Particle Physics. In this convention, the 7r + 7r + scattering length is positive, corresponding to a repulsive 
interaction. 



TABLE I: The parameters of the MILC gauge configurations and domain-wall quark propagators 
used in this work. The subscript I denotes light quark (up and down), and s denotes the strange 
quark. The superscript dwf denotes the bare-quark mass for the domain-wall fermion propagator 
calculation. The last column is the number of configurations times the number of sources per con- 
figuration. For the ensembles labeled with "P±A", propagators that were periodic in the temporal 
direction were computed in addition to those with anti-periodic temporal boundary conditions 
(bm res was not computed for these sets but is expected to be very close to the result from the 
corresponding set of propagators with anti-periodic temporal boubdary conditions). 



Ensemble 


bmi 


bm s 


i dwf 

bm^ J 


i dwf 

bm s J 


10 3 x bm res 


# of propagators 


2064f21b676m007m050 


0.007 


0.050 


0.0081 


0.081 


1.604 


1038 x 24 


2064f21b676m010m050 


0.010 


0.050 


0.0138 


0.081 


1.552 


768 x 24 


2064f21b679m020m050 


0.020 


0.050 


0.0313 


0.081 


1.239 


486 x 24 


2064f21b681m030m050 


0.030 


0.050 


0.0478 


0.081 


0.982 


564 x 20 


2896f2b709m0062m031 


0.0062 


0.031 


0.0080 


0.0423 


0.380 


1001 x 7 


2896f2b709m0062m031 P±A 
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1.552 


137 x 8 


2864f2b676m010m050 P±A 
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0.0138 


0.081 
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274 x (2+2) 



IV. METHODOLOGY AND DETAILS OF THE LATTICE CALCULATION 

The computation in this paper uses the mixed-action lattice QCD scheme developed by 
LHPC [TBI HE] and described fully in Ref. j6] where multi-pion systems are investigated in 
detail. The present calculations were performed predominantly with the coarse, asqtad- 
improved [2Ql El] MILC configurations generated with rooted staggered sea quarks [22J 5 
with a lattice spacing of b ~ 0.125 fm, and a spatial extent of L ~ 2.5 fm. These con- 
figurations were HYP-smeared [231 EU E3 |26] and used to generate domain-wall fermion 
propagators [27J ESI EH1 EH1 |3T] . The strange quark was held fixed near its physical value 
while the degenerate light quarks were varied over a range of masses 6 ; see Table [I] and 
Refs. [31 E21 E31 EH ES] for details. Further calculations were performed using one fine MILC 
ensemble, with b ~ 0.090 fm and L ~ 2.5 fm and on a larger volume coarse MILC ensemble 
with b ~ 0.125 fm and L ~ 3.5 fm. On the coarse, 2.5 fm MILC lattices, Dirichlet boundary 
conditions were implemented to reduce the original time extent of 64 time-slices down to 32 
and thus save in computational time. While this procedure leads to minimal degradation 
of baryon signals, it does limit the number of time slices available for fitting meson proper- 
ties. For the fine MILC ensemble and the large volume coarse MILC ensemble, anti-periodic 
boundary conditions in time were implemented and all time slices are available for analysis. 
In order to extract precise measurements from the configurations, propagators were com- 
puted from multiple sources displaced both temporally and spatially on each configuration 
The correlators were blocked so that one averaged correlator per configuration was used 
in the subsequent statistical analysis. Further blocking over sets of configurations did not 



5 The results of this paper assume that the fourth-root trick recovers the correct continuum limit of QCD. 

6 The mass of the valence quark was chosen to yield a pion mass that matches the pseudo-Goldstone 
staggered pion [6J. 



change the resulting uncertainties associated with each observable. 

For the fine MILC ensemble and the large volume coarse MILC ensemble, additional 
propagators were computed using periodic temporal boundary conditions for a subset of the 
source points. Combining periodic and anti-periodic (in time) propagators from the same 
source leads to a propagator that is periodic in time over twice the temporal extent of the 
lattice. These ensembles, labeled "P±A" in Table [TJ served to explore the effects of temporal 
boundary conditions. 7 

In order to determine the interaction energy in multi-meson systems, both the single- 
meson, Ci(t), and n-meson, C n (t), correlation functions were computed, where t is the 
Euclidean time separation between the hadronic source and sink operators. For simplicity, 
meson correlation functions involving only zero-momentum states were calculated. The 
single-kaon correlation function is 

C K -(t) = J2(K + (t,x)K~(0,0)}, (17) 

X 

where the sum over all spatial sites projects onto the zero-momentum state, p = 0. A 
correlation function which projects onto the n K~ ground-state is 



C nK -(t) = (( y £ i K+{t,x)\ ( K-(0,0)) ) 



(18) 



which can be constructed from one light and one strange quark propagator. In eqs. (17) and 



(18), K + (t,x) = s(t,x)7 5 it(t, x) is a Gaussian-smeared interpolating field for the charged 
kaon. In the relatively-large spatial volumes used in the calculation, the interaction energy is 
a small fraction of the total energy that is dominated by the mass of the kaon. To determine 
this energy, the ratio of correlation functions, 

GnK- (t) EE ^*- (t) — A e" A ^ 4 , (19) 

[ ^K- {t) \ 

was constructed, with the arrow denoting the large-time, infinite-number-of-gauge- 
configurations limit (far from the boundary). Due to the (anti-)periodic boundary-conditions 
imposed on the propagators computed on the fine, and large- volume lattices, the K~ correla- 
tion function becomes a single cosh function far from the source, while the n-kaon correlation 
functions become sums of multiple cosh's, leading to a non-trivial form for G n K-(t); further 
details are given in Appendix |Aj As an alternative method to calculating the interaction 
energy (and a check of the systematics), a Jackknife analysis of the difference between the 
energies extracted from the long-time behavior of the multiple- and single-kaon correlation 
functions individually was performed, finding results in agreement with those determined 



from eq. (19). The interaction energy is related to the n-particle energy eigenvalues, E n x, 
and the kaon mass, 

AE nK = E nK - n m K . (20) 



7 For the correlators of systems containing multiple hadrons (which are extremely sensitive to truncation 
errors in arithmetic operations), we find that performing the P±A additions and subtractions in machine 
(double) precision leads to mild degradation of signals that worsens with increasing n. This indicates 
that differences in the iterative (conjugate gradient) propagator inversion for the two types of boundary 
conditions are becoming significant. 



A comprehensive review of the lattice techniques we have used in this work can be found in 
Refs. [HI [36]. 
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FIG. 1: The energy-differences (in lattice units), b AE n x, for the multi-kaon systems with 
m n ~ 291 MeV. The lighter (darker) shaded regions correspond to the statistical and systematic 
uncertainties combined in quadrature (statistical uncertainty) obtained from fitting the plateau 
regions. The green line denotes the central value of the fit. 

To indicate the quality of the results of this calculation, in fig. [T] we show the effective 
energy differences, AE n x, as defined in eq. (20), for the coarse 2.5 fm ensemble with m^ ~ 



291 MeV. For each n there is clear, multi-time-slice, plateau from which to extract AE n x- 
Figure [T] shows that the systematic uncertainties in the plateau region tend to become more 
significant as the number of kaons in the system is increased, likely due to the fact that 
the same propagators are being raised to high powers. In the extraction of AE n x from the 
correlation functions, a correlated x 2_ fit is performed over the fitting interval. Further, in 
the extractions of the two- and three- body scattering parameters, correlations among the 
different AE n x are included by performing a coupled n- and t- correlated fit. Systematic 
uncertainties are assigned to these extractions by varying the ends of the fitting intervals. 



8 



It is interesting to note that the uncertainty in the effective n-kaon energy grows with time. 
This is expected when m^ + m^ < 2m K (as is the case for the ensembles studied here) based 
on the general arguments given in Ref. [HT 



V. THREE-MESON INTERACTIONS 



The three-meson interactions extracted from our calculations, m n f* r] 37T - 



are shown in Fig. 2[ in units of the NDA estimate. 
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l/(mxfjc)- Table O contains the extracted values of the two-pion and two-kaon scattering 
parameters and the three-pion and three-kaon interactions for each ensemble. The results of 



1.5 



l.O 



lie- n - 



0.0 



-0.5 L 
1.0 





1.0 




_ 




0.5 








- 






M 




II 


- 


lis? 


0.0 










II 


1 


1 


1 1 


1 




" 




-0.5 








- 




-1.0 


- 



1.5 2.0 2.5 3.0 

m„lf n 



3.5 4.0 



3.4 3.6 
m K /fx 



3.8 



4.0 



FIG. 2: The three-body interactions, m n f^ r/ 37r - and mx/x V3K~ > defined in eq. ( 10 ) , as a function 
of run/ fn and mx/Ix, respectively. The statistical uncertainties and the statistical and systematic 
uncertainties combined in quadrature are shown as the inner and outer error-bounds, respectively. 
The dashed vertical lines indicates the physical point. 



the calculation of the three-pion interaction on the coarse MILC lattices have been presented 
previously, and the result from the fine ensemble is consistent with the closest mass coarse 
results, as can be seen in fig. [2} The three-kaon interaction is found to be consistent with 
zero for all the kaon masses that have been explored, but the uncertainties are moderately 
large. This result is consistent with the trend found in multi-pion systems where there is 
a hint that the three-pion interaction is decreasing with increasing pion mass. As is the 
case for the three-pion interaction [6], it is not yet possible to extract the RGI three-body 
interaction, fj 3K or the underlying parameter t]zk{ia) because of the current uncertainties 
in the calculation. The uncertainties of the three-body interactions measured on the large 
volume ensemble are significantly larger than those shown in the figure and we do not display 
them. 



TABLE II: Extracted values of the two-body and three-body interactions. Note that r] 37T - and 
rj 3K - depend on L, eq. (13) and eq. (14) and are quoted at their physical volumes. The systematic 



uncertainties in the pionic quantities from the coarse L ~ 2.5 fm ensembles have been symmetrized 
about the central value [6]. The pion and kaon scattering lengths obtained on the coarse L ~ 2.5 fm 
ensembles are in agreement with those of Ref. [38, 39 , but have slightly smaller uncertainties 
stemming from increased statistics and from the n-correlated analysis performed here. Quantities 
extracted from the correlators formed from the P±A propagators are consistent with those from 
the individual propagators. 
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VI. THE EQUATION OF STATE AND CHEMICAL POTENTIALS 



The energy of the n K system as a function of volume and of the number of K 's is given 



explicitly in eq. (10) in the large- volume expansion. From the equation of state, the K 



chemical potential is 



Hk- 



dE, 



nK 



dn 



(21) 



y=const 



and can be constructed analytically from eq. (10) or numerically from the results of the 



lattice calculations by using a simple finite difference approximation. The resulting ratios of 
the iT~-chemical potentials to the kaon mass on each of the coarse ensembles are shown in 
Fig. K^as a function of the i^~-density, p%- in units of (2.5 fm) 3 . For comparison Fig. Sshows 
the isospin chemical potential computed from the n tt + systems, expanding upon the results 
of Ref. [B] by including the large volume coarse ensemble. Finally, Fig. [5] shows the kaon and 
pion chemical potentials for the fine ensemble, and a comparison between results obtained 
with anti-periodic boundary conditions and those computed with the P±A propagators. 
We note that for the 12-K~ system, the number density is p K _ = 12/L 3 = 0.77 fm" on 
the ensembles with L ~ 2.5 fm. The solid curves in Figs. [3j [4] and [5] correspond to the 
prediction at 0(L~ 7 ) from eq. (10) using the extracted two-body and three-body scattering 
parameters, and differ insignificantly from that at 0(L~ 6 ). The dotted curves result from 
setting the three-body interactions to zero. 

The results of the lattice QCD calculation are consistent within uncertainties with the 
prediction of tree-level SU (3) l<S> SU (3) r x?T which are shown as the dashed lines in Figs. |3j 
|l]and[5J Given that one expects the tree-level result to be accurate at the ~ 25%-level, this 
is quite a remarkable result. The contributions from the higher-order counterterms and loop 
diagrams have conspired to produce a shift from tree-level that is much less than one expects 
from NDA. 

In comparing these lattice QCD calculations of the chemical potential to those calculated 
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FIG. 3: The K~ chemical potential, //r-- , as a function of the density of K~ 's on the coarse MILC 
lattices. The rectangular regions are the results of the lattice calculation where a finite-difference 
has been used to construct the derivative with respect to the number of K~ J s. The lighter (green) 
rectangles are obtained from the lattices with spatial extent L ~ 2.5 fm, while the darker (red) 
rectangles are obtained from the lattices with L ~ 3.5 fm. The dashed curves correspond to tree- 
level ;\PT. The darker solid curve corresponds to the analytic expression for the energy of the 
ground state in the large volume expansion, eq. (10), using the fit values for a K ~^ and rj- iK -. The 



lighter dotted curve shows the analytic forms with the fit value of a KK but with r/ 3 ^-set to zero. 



in xPT in Section [TTJ we note that they are performed in finite volumes and at non- vanishing 
lattice spacings. These approximations modify the chemical potential expected in the lattice 
calculation from chiral perturbation theory at most at NLO in the chiral expansion. The 
fact that LO xPT is in good agreement with our numerical results, on both the larger 
and smaller volume coarse ensembles and the fine ensemble, suggests that finite volume 
and lattice spacing effects are small. For m n ~ 352 MeV, the two calculations at different 
volumes agree within their uncertainties, also supporting this claim. 

VII. CONCLUSIONS 

Kaon condensation may play an important role in the evolution of supernovae pQ. The 
theoretical analysis of the condensation mechanism presently relies upon chiral perturbation 
theory to determine both the interactions of kaons with the baryonic matter and the kaon 
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FIG. 4: The isospin-chemical potential, ///, as a function of the density of 7r~'s on the coarse MILC 
lattices. The rectangular regions are the results of the lattice calculation where a finite-difference 
has been used to construct the derivative with respect to the number of 7r~'s. The lighter (green) 
rectangles are obtained from the lattices with spatial extent L ~ 2.5 fm, while the darker (red) 
rectangles are obtained from the lattices with L ~ 3.5 fm. The dashed curves correspond to tree- 
level ;\PT. The darker solid curve corresponds to the analytic expression for the energy of the 
ground state in the large volume expansion, eq. (10), using the fit values for a^^ and r/ 37r -. The 



lighter dotted curve shows the analytic forms with the fit value of a\m but with ry 37r -set to zero. 



self-interactions. We have performed the first Lattice QCD investigation of the charged kaon 
condensate, isolating part of the mechanism that may be in action in the interior of neutron 
stars. In particular, we have explored the K~ chemical potential as a function of density. 
We have also updated our previous analysis of the charged pion condensate. Surprisingly, 
LO xPT is found to reproduce the results of the Lattice QCD calculations to within the 
uncertainties of the calculation. This is significantly better than the ~ 25%-level expected 
from dimensional analysis. However, this is consistent with the findings of other lattice 
calculations that higher orders in chiral perturbation theory make contributions that are 
much smaller than naively expected, e.g. Ref. [36], in this range of meson masses. 

To move toward exploring the possibility of kaon condensation in a baryonic background, 
lattice calculations of systems containing multiple nucleons and multiple kaons will be re- 
quired. Only one dynamical lattice QCD calculation of the two-nucleon system has been 
performed to date [32] , and exploration of the multi-nucleon-multi-kaon systems are left for 
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FIG. 5: The K~ chemical potential, (ix-i StS a function of the density of K _, s on the fine MILC 
lattice ensemble (left panel) and the isospin chemical potential, [ii as a function of the density of 
7r + 's on the same ensemble (right panel). The rectangular regions show the results of the lattice 
calculation where a finite-difference has been used to construct the derivative with respect to the 
density of mesons. The lighter (blue) rectangles are obtained from purely anti-periodic temporal 
boundary conditions on the quark propagators, while the darker (purple) rectangles are obtained 
from the P±A propagators. The dashed curves correspond to LO xPT. The darker solid curve 
corresponds to the analytic expressions for the energy of the ground state in the large volume 



expansion, eq. (10), using the fit values of the two-body and three-body interactions. The lighter 



dotted curve shows the analytic forms with the fit value of the two-body interaction but with 
three-body interactions set to zero. 

the future. In addition to the increased exponential degradation of the signal-to-noise ra- 
tio associated with baryonic correlation functions [37J, "disconnected diagrams" (in which a 
quark is created and annihilated at the source) will contribute to the multi-nucleon-multi-fC" 
correlation functions, and therefore the computational resources required to perform such 
calculations will be considerably larger than currently available for the study of multi-hadron 
systems. 
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APPENDIX A: EFFECTS OF TEMPORAL BOUNDARIES 

For temporally periodic gauge configurations and (anti-) periodic quark propagators (or 
the P±A propagators), we are in a position to account for the effects of the temporal 
boundary condition. As is well known, a single meson correlation function, generated from 
propagators that are subject to periodic or anti-periodic boundary conditions, does not 
decay exponentially but is periodic in time and is given by the sum of forward and backward 
propagating exponentials, leading to the behavior 

d(t) ^i A 1 [e-^ + e -^ T -*)] , (Al) 

where T is the temporal period and the source is defined to be at t = 0. For a two 
meson correlation function, there are two contributions; the two mesons can propagate 
either forward or backward or one meson can propagate forward and one backward. This 
leads to a correlator 

C 2 (t) ^i A 2 [e- E * + e -^ T -')] + B 1A , (A2) 

where the second term is t independent as it is a product of exp[— E\t\ and exp[— E\(T — 1)\. 
This can further be extended to the n meson case, where the n-meson correlation function 
is given by 

LfJ 
Cn (t) *»i A n [e- E - 1 + e -^ T -')] + V B n , rn [ e -E n - m t e -E m (T-t) + e -E m t e -E n _ m (T-t^ 



ra=l 



(A3) 



We have used the fact that a positive definite transfer matrix guarantees that the energies 
occurring in these expressions do not depend on the correlator they appear in, although the 
normalizations Bij do. This is only true over physical length scales for the DW fermion 
discretization used herein. 

In principle, these expressions can be used to extract the energies, Ei, . . . , E 12 from the 
measured correlations. However, there are increasing numbers of parameters to determine 
as n increases (e.g., n — 12 involves eight new parameters) and we have only been able to 
successfully perform this analysis up to n = 6. The energies determined with this procedure 
are consistent with those extracted from the plateaus in the effective mass plots of the 
individual correlation functions, but with somewhat larger uncertainties. 

As a further check of our primary analysis, the correlators generated by combining the 
propagators subject to periodic boundary conditions with those subject to anti-periodic 
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boundary conditions, denoted by P±A, (whose behaviors are described by eqs. (Al), (A2) 



and (A3) in which T is replaced by 2T, significantly reducing the contributions of backward 
propagating states) were compared with the corresponding correlation functions using only 
anti-periodic boundary conditions for the quark propagators. On both the fine and large 
volume coarse MILC ensembles, statistical agreement was found for the extracted n-ir and 
n-K energies for all n studied. Fig. [5] compares the two extractions of the n~ and K~ 
chemical potentials for the fine ensemble. 

In the case of a Dirichlet boundary in time, as employed in the calculation of quark 
propagators on the coarse L ~ 2.5 fm ensembles, the nature of the states reflected from the 
walls is less clear and we do not attempt to account for them, instead restricting our analysis 
to a region that is well separated from the walls. 
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